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Abstract 



We apply stochastic quantization method to the bosonic part of IIB matrix model, 
i.e., a naive zero volume limit of large N Yang-Mills theory, to construct a collective field 



> 

X 

■ theory of Wilson loops. The Langevin equation for Wilson loops can be interpreted as 
the time evolution of closed string fields. The corresponding Fokker-Planck hamiltonian 
deduces a closed string field theory which describes interacting Wilson loops with manifest 
Lorentz invariance. 
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It has been conjectured that non-perturbative definitions of superstring theories are 
given by supersymmetric matrix models]l||[^P|. The proposal is strongly supported by 
the fact that the IIB matrix model almost recovers the IIB light-cone string field 
theory^ at the scaling limit 0. Although the light-cone setting is preferable to prove 
the equivalence, we hope to take the continuum limit of IIB matrix model preserving the 
manifest Lorentz invariance for studying its dynamical behaviour such as the spontaneous 
break down of the Lorentz invariance in the context of the superstring field theory. One 
way to keep the manifest Lorentz invariance is to apply stochastic quantization method 
(SQM)0 to IIB matrix model. In this note, we consider only the bosonic part of IIB 
matrix model, which was once discussed to map the large N gauge theories to a string 
, to derive a hamiltonian which describes the time development of Wilson loops and 
to clarify how to take the continuum limit of it. The dynamical aspects of bosonic IIB 
matirx model have been studied]^ [|1^] especially in relation to the U{lY symmetry in the 



Eguchi-Kawai model |Tl|. It is pointed out that the bosonic IIB matrix model describes 
the weak coupling region in the original Eguchi-Kawai model where the U{lY symmetry 
is spontaneously broken p!^, however, the eigenvalue distribution does not collapse to a 



point but has a finite extent at the t'Hooft limit 0. As a naive zero volume limit of 
SU(N) Yang-Mills theory, the partition function of the bosonic IIB matrix model is finite 



and the model appears to be well-defined in a mathematical sense |13|. Our strategy 



in constructing the string field theory is similar to that applied to the minimal matter 



matrix models ||T^ |T^. We mainly use SQM to construct a collective field theory of a 
gauge invariant observable, Wilson loops, in bosonic IIB matrix model so that we do not 
need any gauge fixing procedures. The Fokker-Planck hamiltonian defines a string field 
theoretical description of Wilson loops at the scaling limit. 

We begin with the action of the bosonic part of IIB matrix model, 

S=-^M[A„A^r), (1) 

where A^'s are hermitian N x N matrices. We define the time evolution of the fun- 
damental matrix variables, A^{t + At) = A^{t) + AA^{t), in terms of Ito's calculus 



with respect to the discretized stochastic time r. The time evolution, which recovers the 
probabihty distribution e~'^ with S in the equation (1) at the infinite stochastic time 
hmit r ^ cxD, is given by the following basic Langevin equation, 

A^.(r) = -^[A., [A„ A1](r)Ar + A^,(r) . (2) 

Here the white noises, A^^'s, are also hermitian N x N matrices. Their correlations are 
defined by 

< A^f,{r)ijA^^{T)ki >5= 2ATS^^5iiSjk . (3) 

The Wilson loop is a gauge invariant observable which is calculable in SQM without 
gauge fixing procedure, 

^ M 

Wuih, kM) = -Tr n Un , (4) 

n=l 



where Un = exp{iek^Af^/\' N) . It satisfies an equation, 

KK^M + 0{e') = , (5) 

where A;^ = likl^ + k^+i) , and Vj] = d / d{iek!^) — d / d^iek^+i) . As e — > , we will consider 
the continuum limit such as {ki, k!^, k^j) — > A;^(cr) , and 



n 



d{itk'f,) 6kt'{a) ' 
V^, = X„\i-X„^^eX"'(a) , (6) 

by introducing the world sheet coodinate a, — tto; < cr < vro; , with the "length "parameter 
of the Wilson loop, a = Mejlv: . Then the equation (|) corresponds to, k^{a)X' ^^(oyWoSkio^ 
, which ensures the reparametrization invariance of the Wilson loop with respect to the 
world sheet coordinate (t[§]. 

The generalized Langevin equation for the Wilson loop is derived from the variation 

1 M 

AWu = - E Tr(f/i...Af/„...f/M) 

m=l 



+ ^ E Tr(f/i...At/„...Af/„...f/M) + 0(Ar3/') . (7) 



l<m<n<M 



From the basic Langevin equation (0) and the noise correlation (|^), we have, 

Jo 

- 2e2Arl(A;C)2 [' ds f AUn{l - t)Trf/„(t) + 0{At^'^) , (8) 
iV JO Jo 



where Un{s) = exp(isY'„) , and Yn = ekl^A^/ y N . The first term in r.h.s. of (0) is 
equivalent to the insertion of the equations of motion and the noise variables which can 
be evaluated by the similar procedure in Ref. 0. Then we obtain 

1 A/ 1 M 

-J2MUi...AUn...UM) = -Y.^T{U^...Un{^eKA^,/vT 

n=l n=l 

- Ar^(Vi:)^-Are2(A:^)%}[/„+,...[/M) , (9) 

where In = jj Jods J^dUn{l — t)Trf/„(t) . The second term in r.h.s. of (0) will give a 
splitting interaction of the Wilson loop. It reads, 

^ J2 Tr([/i...Af/„...Af/„...f/M) + 0(Ar3/2) 

-'\<m<n<M 

= -2Are2^ ^™ " A;„Tr([/i...f/^f/„+i...[/M)Tr(f/^+i...f/„) . (10) 

l<m<ra<M 

We here notice that the splitting interaction terms make effects to renormalize the ma- 
trices, (/„)jj's, in equation (j^) at the continuum limit. A possible contribution may come 
from, X]i<n<Af Tr(f/i...Af/„Af/„_i_i...f/A/) , and yields terms which renormalize the kinetic 
term k^{aY at the continuum limit without changing the sign of it. In the following, we 
assume that {In)ij = I^ij and the factor / includes these contributions. 

Now we take the continuum limit. We introduce the renormalized (continuum) 
stochastic time, dr = e^Ar , and the renormalized Wilson loop 

1 

WM^)] = ^~'Wm = ^Tr n Un . (11) 

n=l 



To take the well-defined continuum limit e ^ in the generalized Langevin equation 
for Wilson loops, we require that the free part except the noise term in equation (^) is 
identified to the constraint of the reparametrization invariance with respect to the time 
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coordinate r ; (j€^^{a)^ + ■ We also require for the splitting interaction 

part, equation ([To|), to be finite at the continuum limit. 
From these requirements, we assume 

gle""-^ = a'\le^-'' = 1 ,6 = 1 + D , (12) 

up to some numerical constants. Then we have, 

- dr /_™ (^X'^(a))'|l^„[A;(a)] 

where a = ai + a2- The overlap condition for the splitting interaction is the following, 

A;(3)(a3) - 0ikM\ai) - 92k^^\a2) = 0i = 0{7iai - \a\) 

02 = 9i\a\ - 7ra,) 
6; = 01 + 02 = I 

—Tcai < cr < Txai 
nai < a < 7r(ai + 02) 
(7 + Trai — 7r(ai + 0:2) < o" < —ttoi 

as = a + TT/Ss -7r(ai + 0:2) < o" < TT{ai + 0^2) 

- "2 < /33 < "2 (14) 

The parameter is introduced to represent the splitting point which is specified by 
0-3 = ±nai + on the Wilson loop. We denote the overlap condition as 5(1, 2, 8^3). 
The derivative coupling at the splitting interaction point is defined by kj^\±'7iai + TrP3) = 
^{k^\±7rai + 7cP3 + e) + k'^^\±7cai + 7cP3 — e)} . Since this quantity diverges as e~^/^ at the 
interaction point 0, the factor e^^/^ is necessary in the equation (|T3|) for the well-defined 
continuum limit. 
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The noise variable is also defined by, 

1 M 

A^W^a)] = -— ^ IV (c/i • • • [/„ {iek>;,A^jVN} U^+i ■ ■ ■ Um) (15) 

n=l 

We obtain the following correlation of them at the continuum limit. 



= -2dTG.. 



St 



-Ql J —0.2 J 



x{6i/2fc«(7rai+7rA)}{e^/2^(2)(7r/52)}5(l/3„2^„3)iy„3[A;(3)(a3)] (16) 
where a-^ — ax^ai- Here we have identified 

We notice that, for the merging interaction, the overlap condition is slightly different 
from that of splitting interaction. 

A;(3)(a3) - ^lA;^^^!) - ^2A;('n^2) = 

(Ti = cr + TT/Si — Trai < cr < Trai 

cr + 7r/32 — Trai Trcti < (7 < 'K{pL\ + q;2) 

cr + 7r/32 + Trai — 7r(Q;i + a-i) <o< — Trcci 

(73 = (7 — 7r(Q;i + a'l) <o< 7r(ai + a-i) 

— OL\<i5\< ai, - Q;2 < /?2 < Q;2 (18) 

The overlap (5-function is now denoted by 5(1/3^,2^2,3) , The merging interaction point 
is specified by o"i = ±7rQ;i + 7r/3i and (T2 = 7r/52 • The derivative coupling at the merging 
interaction point is also different from that at the splitting interaction point. They are 
given by 

A;W(7rQ;i + 7r/?i) = ^(A;(3)(7rQ;i - e) + ^(^^(-Trai + e)) , 

^(2)(7r/32) = ^(A;(=^)(7rQ;i + e) + V<^\-T^a^ - e)) . (19) 
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Before deriving the Fokker-Planck hamiltonian, we add some remarks on the scahng 
hmit which is defined by equations ( |T^ ) and (p!?]). Since we fix the string tension a' and 



the string couphng constant Gst to be finite at the large N hmit, the equations ([T2|) and 
(p!7D reads, 

gl ^ N-^^-D)/{i+D) (20) 



A weU-defined continuum hmit requires e ^ at ^ oo which ensures the scahng 
dimension of the stochastic time r to satisfy 1 + D > . We here require a stronger 
condition D > to obtain the smooth continuum hmit for the time development of 
the stochastic process. By redefining the fundamental matrix variable to Afj_\/N 
both in the action (||) and the definition of the Wilson loop (^, the coupling in 
(|1]) is replaced by g^/N which is equivalent to the original coupling constant of JIB 
matrix model g^. According to the correspondence g^ = g'^N, the space-time extent is 
expected to be characterized by the quantity R ~ ^N^^^ = y/g^ p. It seems that, if a 
continuum limit may exist, there are three possible cases, 0<-D<l,D = l and D > 1, 
which distinguish the extent of eigen value distribution, i.e., the space-time extent. For 
< D < 1, since R at A^ — > oo , the target space-time is collapsed. The case 
D = 1 in (|1^) realizes the t'Hooft limit and it deduces A^e^ = 1/Gst- The condition 
was discussed in the light-cone setting]^. In this case, the space-time may have a finite 
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extent even at the large N limit. Then we may consider a dual picture to interpret the 
momentum k^{ij) for the closed string by assuming the boundary condition for the target 
space time coordinte [0, n^R = x^{27r) — a;^(0) , with n^, the winding number in the n 
direction. The last possibility is to consider the limit R —>■ oo which means D > 1. In any 
cases, a remarkable fact is that since Ne^ = l/y/G^t is fixed, the intrinsic Wilson loop, 
TrH^lif^n 5 is not scaled at the scaling limit. Our analysis is restricted to the bosonic 
part of IIB matirx model, however, the scaling relations, ([T2| ) and (p!7D, hold even for full 
IIB matirx model. The dynamical behaviour of the bosonic IIB matrix model is expected 
to be quite different from that of full IIB matrix model and it should be determined from 
the large N behaviour of the numerical factor / e^^^. If we assume that the factor I 
is a function of e,gQ,N] I = I{e, g^, N), then the condition to determine the index D is 
given by gQl{N^^^^^^^\ Qq, N) ^ 0(1) . Further investigation of the dynamical behaviour 
of this numerical factor, however, is beyond the scope of this note. 

In our approach, the stochastic process of the Wilson loop, which is described by 
the generalized Langevin equation, is interpreted as the time evolution in a string field 
theory. The corresponding Fokker-Planck hamiltonian defines a bosonic closed string 
field theory. In terms of the expectation value of an observable 0{Wa), a function of 



WaS, the Fokker-Planck hamiltonian operator Hpp is given by[|U 

<W^{0)\e-^^^''O{W^)\0>=<O{Wa{r)^)>^ . (21) 

In r.h.s., Wa{T)^ denotes the solutions of the Langevin equation with the appropriate 
initial configuration Wa{0). In l.h.s., Hpp is given by the differential operator in the 
well-known Fokker-Planck equation for the expectation value of the observable 0{Wa) 
by replacing the closed string variable lVQ,[/i;((j)] to the creation operator 14^q,[A;((j)] and 
the differential 5/5WQ,[/i;((T)] to the annihilation operator nQ,[A;(cr)]. The commutation 
relations are 

[tla[kia)],W^,[k\a')]] = Sia - a') [] m^)-k'ia)) . (22) 

—7va<a<7Ta. 
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At the scaling limit, we obtain the continuum Fokker-Planck hamiltonian, Hfp-, 
Hfp = da jDk{a) | (^^X'^(a))'| l^„[fc(cr)]n4A;(a)] 



+ 



/"OO /"ft /'CV2 P 

2 da dai d/33 pA;W(ai)m(2)(cr2)m(3)(a3){ei/2A;(3)(7rai + 7r/33)} 

Jo Jo J-a2 J 

roo roo pc^i POi2 P ^ 

+ Gst / dai / da2 / df3i / rf/?2 bk^^\ai)Vk^^\a2)Vk^^\a3){e^/^k^^\TTai + 7r(3,)} 

Jo Jo J-a-i J-a2 J 

■{ei/2^(2)(7r/52)}5(l/3„2^„3)iy„3[A;(=')(a3)]n„Jfc«(ai)]n,JA;(2)(a2)] (23) 

By redefining the momentum integration variable, k^^\as) = k^^^as) = k^^\a + TxjSz) 
with 0-3 = 0", (— 7r(ai + 0:2) < o- < 7r(ai + 0:2)) , we can eliminate the /53 dependence in the 
splitting overlap, 5{l,2,?)p^) , and in the derivative coupling, A;*^^)(±7rai + vr/53) . We can 
also eliminate the /5i and /?2 dependence in the merging interaction. After the integration 
with respect to the angle variables j3i, {i = 1, 2, 3), the Fokker-Planck hamiltonian finally 
takes the form, 

,2/1 



^FP = |^~f/a^A;(a)|™ rfa|(A;'^(a))% (^^X'^(a)) |H^„[A;(a)]n^ 

+ 2 Tda r dai{2Tia2) fDk'^^\a^)Vk^^\a2)Vk^^\a:i){t^/^U^\Tiai)] 
Jo Jo J 

.{eV2^(3)(_^^^)|5(1^2,3)iy„JA;«(ai)]iy„Jfc(2)(a2)]n,3[^^'H^3)] 

+ Gst /"c^ai /"c/«2(47r2«i«2) fDk^^\ai)Vk^^\a2)Vk^^\a3){e'^^k^^\na,)} 
Jo Jo J 

■{ei/2^(^)(0)}5(l,2,3)iy«3[A;(3)(a3)]n,JA:«(ai)]n,JA;(2)(a2)] (24) 



where the overlap 5-function 5(1, 2, 3) is given by the equation (|T4D with = 0. Remark- 
ably, the hamltonian consists of the splitting and the merging interaction terms where 
only three strings interact. At the weak string coupling limit Ggt —>■ 0, which is equivalent 
to the large N limit with finite e, the merging interaction vanishes. 

In the definition of the string field theory interpretation (|21|), the initial state is defined 
by < IVq(0)| =< 0\exp{ J^daVk{a)Wa[k{a)]fla[k{a)]} , where the initial configuration 
must be the solution of the large N factorized Schwinger-Dyson equation ( or the loop 
equation )|1^. The S-D equation is given by the limit AWa — > in the generalized 
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Langevin equation (jT3|), 

+ 2 f'dai r'dp, fofc<"(ffi)I'*:'"(<^2){6"'i''"(irai + ir/33)} 

JO J-a2 J 

.(eV2p)(_^tti + 7r/33)}5(l,2,3ft) < PF„Jfc«(ai)]W^„JA;(2)(a2)] >= . (25) 

It is clear from the correlation in ([T6D that the leading finite N effect on the correlation 
functions is given by the order of 1/A^^[0 and the S-D equation is factorizable at large 
N. 

In conclusion, we have derived the continuum Fokker-Planck hamiltonian from the 
stochastic process defined by the bosonic part of IIB matrix model. The structure of 
the splitting and the merging interactions we have derived are covariantized versions of 
those in the bosonic part of the type IIB Green-Schwarz superstring field theory in the 
light-cone gauge. The Fokker-Planck hamiltonian may be regarded as a closed string 
field theory hamiltonian which describes the time development of Wilson loops. We have 
constructed the field theory hamiltonian such that it does not have a free theory and the 
splitting interaction remains in the weak string coupling limit Gst . The situation is 
analogous to the string field theories constructed from the old-fashioned matrix models. 

Now we comment on some issues which remains to be studied. If we assume that 
the bosonic IIB matrix model describes the weak coupling region where U{lY symmetry 
is broken spontaneously, the Schwinger- Dyson equations in this model do not recover 
the planar perturbation theory defined by the continuum Yang- Mills theory]!^. We also 



do not have evidence that implies the equivalence of the Fokker-Planck hamiltonian to 
collective field theories of Wilson loops in the canonical formalism of standard lattice 
gauge theories [p!^] or to the bosonic closed string field theory In a formal sense, the 



Fokker-Planck hamiltonian is a manifestly Lorentz invariant hamiltonian constraint for 
observables such as correlation functions of Wilson loops at the infinite stochastic time 
r ^ oo. 

<i/^p{iy,Jfc«(a)],iy,JA;(^)(^)],-,W^a„[A;(")(a)])} >= . (26) 
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Thus, in light-cone variables, the constraint may be reduced to a light-cone hamiltonian 
as expected from the fact that the IIB matrix model almost recovers the light-cone IIB 
superstring field theory [^]. Another possibility is, by reversing the operator ordering of 
the Fokker-Planck hamiltonian, we obtain the Lorentz invariant Fokker-Planck equation 
for the probability distribution functional of the Wilson loop P(Wa[k{a)]). 

^PiW^a)]) = -HMW^iHa)], ^^^r|(^)^(^«[M^)]) • (27) 

The solution of this constraint, i.e., the Boltzmann weight of the path-integral represen- 
tation defined at the infinite stochastic time, may specify a classical action of interacting 
Wilson loops. 

The question whether the present continuum limit is realized dynamically or not is 
yet to be studied, the behavior of the bosonic IIB matirx model may be clarified by 
evaluating the small string contribution to the kinetic term (fc^(cr))^. For IIB matrix 
model, we expect the situation will be changed because of the existence of the stable 
perturbative vacuum and the weak string coupling limit of it may remain a free theory. 
We hope to report on the analysis of full IIB matrix model in future publication. 
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